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The viscous drag on a colloidal particle pulled through solution by an optical trap is large enough 
that on experimentally relavant time scales the mechanical force exerted by the trap is equal and op- 
posite the viscous drag force. The rapid mechanical equilibritation allows the system to be modeled 
using equilibrium theory, where the effects of the energy dissipation {thermodynamic disequilibrium) 
show up only in the coordinate transformations that map the system from the laboratory frame of 
reference, relative to which the particle is moving, to a frame of reference in which the particle is, 
on average, stationary and on which the stochastic dynamics is governed by a canonical equilib- 
rium distribution function. The simple equations in the stationary frame can be analyzed using the 
Onsager-Machlup theory for stochastic systems and provide generalizations of equilibrium and near 
equilibrium concepts such as detailed balance and fluctuation-dissipation relations applicable to a 
wide range of systems including molecular motors, pumps, and other nano-scale machines. 
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It is often emphasized that motion of a small particle 
through water under the influence of an external field is 
a thermodynamically non-equilibrium process. This is 
of course true since there is a net flow of energy from 
the source of the field to the bath in which the particle 
moves and in which the energy is dissipated. Thus it may 
seem that the title of this paper is a nonsequitor - how 
can an equilibrium theory possibly describe a patently 
non-equiibrium process ? 

What is overlooked is that for colloidal particles (radius 
< 10/um) on time scales long compared to the thermal 
(velocity) relaxation time {tv = m/R k, 10~^s) the vis- 
cous drag is so large that acceleration is negligible (the 
over-damped limit) and the viscous drag force is equal 
and opposite the mechanical force P, 0] . This rapid me- 
chanical equilibration allows a very simple description 
of the dynamical behavior of single particles interacting 
with moving potentials. A coordinate transformation is 
used to shift the description of the system from the lab- 
oratory frame of reference, relative to which the colloidal 
particle is moving, to one in which the particle is, on 
average, not moving. The stationary probability density 
function for the position of the particle on this new co- 
ordinate is given by a canonical equilibrium distribution 
function from which all thermodynamic quantities can be 
calculated. Generalization of traditional thermodynamic 
equilibrium and near equilibrium concepts such as de- 
tailed balance and fluctuation-dissipation relations valid 
for systems arbitrarily far from thermodynamic equilib- 
rium, but close to mechanical equilibrium, are also ob- 
tained. 

The specific system discussed in this paper - a colloidal 
particle pulled along with a harmonic well that is caused 
to translate at constant velocity - is a particularly sim- 
ple example that has been studied experimentally j3| and 
theoretically 0, 0| . The general perspective developed is 



that in most cases a nanoscale object in solution under 
the influence of an external force field is itself in equilib- 
rium and undergoes equilibrium fiuctuations in a suitable 
coordinate system. The bath of course is not in equilib- 
rium with the source of the field. The particle serves as a 
"conduit" for energy to fiow from the source to the bath, 
and the energy flow is modulated by the equilibrium fluc- 
tuations of the particle. This perspective is expected to 
apply to a wide range of systems, including not only col- 
loidal particles but also to molecular motors and other 
nano-scale machines 6] . 

Consider a colloidal particle pulled in one-dimension 
through solution by a harmonic potential exerted by an 
optical trap moving at constant velocity (Fig. 1) as in the 
experiment of Wang et al. The equation of motion 
is 



mx 



Rx 



's{x — vt) + e{t) , 



(1) 



where m is the mass of the particle, R is the coefficient 
of viscous drag, s is the spring constant of the harmonic 
optical trap, v is the velocity at which the trap trans- 
lates in space, and e(<) is Gaussian white noise with 
statistical properties < e{t) >= and < e{ti)e{t2) >— 
2RkBTS{t2—ti) where ks is the Boltzmann constant and 
T is the Kelvin temperature. We can remove the explicit 
time dependence with a change of coordinate to 



V{t - Tp) 



(2) 



where Tp = R/s is the positional relaxation time. At 
a = the system is in mechanical equilibrium where the 
force due to the trap is equal and opposite the drag force 
s{x — vt) = Rv when the particle moves at the velocity v 
of the trap. Substituting 



+ V{t~ Tp) 



(3) 
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into Eq. (1) we have 

ma + Ra + sa = e(t) 



(4) 



the equation for a damped harmonic oscillator with noise. 
Subsequent to the velocity relaxation time Ty — m/R the 
acceleration is negligible (i.e., the "overdamped" limit) 
and we have 



Ra + sa — e(t) 



(5) 



Eq. (4) is precisely the equation studied by Onsager and 
Machlup Ig. The conditional probability density that 
the particle moves to position aj at a time t + At given 
that it starts at position ai at time t is 



P{aj,t + At\ai,t) 



exp 



- At/- 



PQi] 



2fcBT(l-e-2A'/T) 



(6) 



Eq. 5 is rather complicated. The ratio of the conditional 
probability densities is much simpler: 



P{aut + At\aj,t) 



exp 



si at 



(7) 



Time {At) has disappeared entirely! Similar rela- 
tions were shown for arbitrary potentials using the On- 
sager /Machlup approach |^. Eq. 6 reflects the fact that 
the system obeys a detailed balance condition on the 
transformed coordinate 



Pcq(ai)-P(aj, i + Ai|Q!i,t) = Poq(aj) P{au t + Aijaj, i) 



The optical potential energy is 
U{x,t) = ^(x - vtf 



2(a 



(8) 



(9) 



from which we see that the term in the exponent in Eq. 
(7) can be written s(a? — a?) = AU{t) + 2RvAa where 
AU{t) = U{xj{t)) - U{xi{t)), and Aa = aj - a-,. Since 
this term depends only on the potential energy at the two 
endpoints, and on the distance between them Aa(= Ax), 
Eq. (7) can be rewritten in the laboratory coordinate as 



P{xj,t + At\xi,t) 
P{xi,t + At\xj,t) 



exp 



-AUit) - 2RvAx 



(10) 



where Aa; = Xj — Xj. The quantity i?u|Ax| is the mini- 
mum work needed to reliably move a distance Ax at con- 
stant velocity v through the viscous medium . Eq. (10) 
is a generalized detailed balance condition for systems in 
mechanical equilibrium but for which no thermodynamic 
equilibrium state exists. For w = Eq. (10) reduces 
to the standard detailed balance condition derived from 
knowledge of the thermodynamic equilibrium to which 
the system relaxes in the long time limit. 



Eqs. (6), (7), and (10) with time broken into incre- 
ments At that are long compared to r„ but short rela- 
tive to Tp can form the basis for a general treatment of 
over-damped stochastic dynamics using the variational 
approach developed by Onsager and Machlup ||7| which 
can be extended to arbitrary potentials [l|, |8j with time 
dependence jld|. 

Here we focus on results that hold for times long com- 
pared to Tp after which the distribution on the trans- 
formed coordinate obeys the canonical stationary equi- 
librium probability density function 



Peq(a) 



^{2^kBT/s) ' 



(11) 



with the mean square displacement from a = given by 
a fluctuation dissipation relation 



{a') = ^s/{2.kBT) / Wexp(-|^) = M: 



and with the two time auto-correlation function [l^ 



:i2) 



ksT 
2s 



exp [-{tl - t2)/Tp] 



(13) 



The distribution function p^q^a) is Gaussian. This is be- 
cause the moving potential is harmonic (i.e., quadratic). 
There is a tendency toward misidentification of a 
Gaussian distribution with a thermodynamically near- 
equilibrium distribution and vice versa. This is not cor- 
rect. For a strongly asymmetric potential such as the 
asymmetric sawtooth often used in discussions of Brow- 
nian motors (lli. .12.. .l.'i] the distribution function is not 
well approximated by a Gaussian even at thermodynamic 
equilibrium. On the other hand, with the moving har- 
monic trap used here, the distribution is Gaussian after 
t > Tp irrespective of the velocity v of the trap, i.e., ir- 
respective of how far the system is from thermodynamic 
equilibrium. The question of whether a distribution func- 
tion for some quantity is or is not Gaussian has as much 
or more to do with the underlying symmetry of the sys- 
tem as it does with whether the system is or is not near 
thermodynamic equilibrium. 

In the experiment of Wang et al. ^ a feedback loop 
was used to assure that the trap moved at constant veloc- 
ity and data were collected only after the average velocity 
of the particle had attained a stationary value (i.e., af- 
ter a time t > Tp). The particle then rides along with 
the moving potential at the average velocity v in the lab- 
oratory frame. The average position in the tranformed 
frame is < a >= with an average power consumption 



Pa 



RV" 



(14) 



The experiment was carried out close to both velocity 
and positional equilibrium. The trap constant used was 
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fairly small (s '-^ 10~^ N/m= 10~^ kg/s^), the radius of 
the particle was r ^ lO^^m so the mass was m ~ Pp ^ 
10~^^ kg (where p ^ lif' kg/m-^ is the density of the 
particle), and the viscous drag coefScient according to 
the Stokes formula was R ~ Q-Krjr ~ 10~^ kg/s, where 
77 ^ 10^"^ kg/(m s) is the viscosity of water. The velocity 
at which the trap was translated was v ^ 10"^ m/s. The 



velocity relaxation time was — m/R 



10 
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the 



positional relaxation time was Tp = R/s ^ lO^^s, and 
the average power input was ^ 2.5fcBT/s. 

As the particle fluctuates about the mechanical equi- 
librium point a = the instantaneous power Pinst needed 
to keep the trap moving at constant velocity fluctuates 



X s{x — vt) = V X s{a — vTp) 



(15) 



where s{x — vt), the reactive force due to the particle 
acting on the potential, is equal and opposite the force 
of the potential acting on the particle. Thus the amount 
of work supplied by the external source in a time At is 

ft+At 

AW = vs I dt'{a - VTp) , (16) 



the mean value of which is 

(AW) = V^STpAt = T'avgAt 



(17) 



There are fluctuations about the mean due to the fluctu- 
ations of the particle. Because a is Gaussian distributed 
W is also Gaussian distributed. Thus we need only to 
calculate the variance 



{AW - {AW))^'j ^vs/ (//+'^* a{t')dt' 



(18) 



= /t*'*''^* /t*'*'^* (a(ti)a(i2)) dtidt2 



It Jt 

to get the probability density function. Using Eq. (13) 
for the correlation function under the double integral in 
Eq. (18), and taking the leading term for large t, we find 
for the variance 



(AW^ - (AT4^))') = 2fcsr7'avgAt , 



(19) 

and the probability density for the work done by the 
source is 



P{AW) 



y/ATTkBTV^vgAt 



(20) 



The work probability density obeys a generalized fluctu- 
ation dissipation relation 



sity 

m 



P{~AW) = P(AVF)exp. 

and the concomitant equality 
/ -AW 



-AW 



(21) 



exp ■ 



kBT 



°° AW 

exp P{AW)dAW = 1 



that follows from the fact that the work probability 
density is normalized J^^P{-AW)dAW = 1. The 
fluctuation-dissipation form of Eq. (19) is obvious from 
the physical perspective - those trajectories in which the 
internal energy of the bath increases are relatively more 
likely than those trajectories in which the internal en- 
ergy of the bath decreases . For over-damped systems 
the generalized fluctuation-dissipation relations follow al- 
most immediately from Onsager and Machlups work 0]. 
Bochkov and Kuzovlev |l5l | used only microscopic re- 
versibility in their derivation, which is consequently more 
general. More recently, equations of the form of Eq. (19) 
have been termed Fluctuation relations llfil and Eq. (20) 
a non-equilibrium work energy relation |l7| . 

From Eq. (19) it is apparent that it is more likely that 
work is done by the source on the bath, than that the 
thermal fluctuations in the bath conspire to do work on 
the source. However, for small At there is an apprecia- 
ble chance that the fluctuations actually reverse the net 
flow of energy and that the work done in the interval is 
negative, i.e., that the bath does work on the source, a 
so-called "violation of the second law" We can cal- 
culate the probability for this to happen by integrating 
the work probability density function from —00 to 0, 

P [AW < 0] = y P{AW)dAW = ierfc ^ 



V At 



(22) 



4fcBr J 

(23) 

A plot of P{AW < 0)/ [1 - P{AW < 0)] against At is 
shown in Fig. 2. The curve agrees well with the ex- 
perimental result of Wang et al. |^. It is not at all 
surprising that such reversals of energy flow occur, but 
it is perhaps counter-intuitive that they are evident for 
micrometer colloidal particles on time scales as long as 
seconds, many times the positional relaxation time of the 
particle. 

By focussing on the mechanically equilibrated motion 
of small scale systems in water simple equilibrium rela- 
tions can be obtained for the probability density func- 
tions in a coordinate system shifted from the laboratory 
frame of reference. These relations contain all informa- 
tion concerning the thermodynamically non-equilibrium 
exchanges of energy between the external source, the sys- 
tem of interest, and the bath. The resulting equations 
provide generalizations of relations derived from equi- 
librium thermodynamics such as detailed balance and 
fluctuation-dissipation relations. The perspective offered 
here is expected to provide important insight into a wide 
variety of nanoscale systems including biologically im- 
portant processes. For proteins we expect for parame- 
ter values s ~ r ^ lO^^m, m ~ lO^^^kg 
and R 10~^^kg/ s from which we find the velocity 
relaxation time t„ ^ 10~^'^s, and the positional relax- 
ation time within a nanometer well, Tp = R/s ^ lO^^s. 
Consequently a theory based on mechanical equilibrium 
is very valuable for understanding the large scale mo- 
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tions involved in protein folding, nanometer stepping of 
biomolecular motors, and large conformational changes 
of signal- and other energy- transducing proteins, all pro- 
cesses that take place while the macromolecule is in me- 
chanical equilibrium but is driven by thermodynamically 
far from equilibrium processes in the environment. 

The generalized fluctuation-dissipation theorems are 
complementary to work on nanoscale brownian motors 
0, 0| and fluctuation driven transport 0] . The gener- 
alized fluctuation-dissipation theorems focus on symme- 
tries that are preserved even in the presence of external 
directed forces that drive the system out of thermody- 
namic equilibrium, while the "brownian motor" mecha- 
nism shows how structural asymmetry can combine with 
directionless energy input and thermal noise to drive di- 
rected motion. A brownian motor can be viewed as a 
molecule or nanoscale device that serves as a conduit 
coupling two external sources to a heat bath in such a 
way that the flow of energy from the stronger source can 
rectify the occasional reversal of the flow of energy be- 
tween the bath and the weaker source, allowing energy to 
be pumped from bath to do work on the weaker source. 
The energy for the reversal is provided by the stronger 
source, but the mechanism takes advantage of the om- 
nipresent fluctuations in the energy flows due to thermal 
or other sources of noise. These fluctuation driven trans- 
port processes |20j share far more in common with the 
coupled transport processes described by Onsager 0, 0| 
than they do with the mechanisms of macroscopic motors 
and pumps. 
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